Quantum discord in nuclear magnetic resonance systems at room temperature 
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We review the theoretical and the experimental aspects regarding the quantification and identifi- 
cation of quantum correlations in liquid-state nuclear magnetic resonance (NMR) systems at room 
temperature. We start by introducing a formal method to obtain the quantum discord and its 
classical counterpart in systems described by a deviation matrix. Next, we apply such a method 
to experimentally demonstrate that the peculiar dynamics, with a sudden change behaviour, of 
quantum discord under decoherence, theoretically predicted only for phase-noise channels, is also 
present even under the effect of a thermal environment. This result shows that such a phenomena 
are much stronger than we could think, at principle. Walking through a different path, we discuss 
an observable witness for the quantumness of correlations in two-qubit systems and present the 
first experimental implementation of such a quantity in a NMR setup. Such a witness could be 
very useful in situations were the knowledge of the nature of correlations (in contrast of how much 
correlations) presented in a given state is enough. 
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I. INTRODUCTION 

A. Bell's Inequalities, Entanglement and the 
No-Local-Broadcasting Theorem 

The importance of the search for a deep understand- 
ing about the quantum nature of correlations among the 
parts constituents of a composite system began to be rec- 
ognized with the seminal works published by Einstein, 
Podolsky and Rosen [l| and Schrodinger [2] in 1935. Re- 
cently, developments in the study of quantum correla- 
tions have been linked with advances in quantum infor- 
mation science. Initially, the only remarkable aspect of 
quantum correlations was their non-local nature, named 
entanglement by Schrodinger and mathematically for- 
malized by Bell in 1964 Q. Bell proved that certain 
relations among the average values of observables of a 
composite system must be satisfied under the assump- 
tions of realism and locality. However, such inequalities 
could be violated for quantum systems in certain entan- 
gled states. These relations are known, in general, as 
Bell's inequalities. States violating a Bell's inequali^ 
are said to possess non-local quantum correlations 

In 1989, Werner ff\ reported a groundbreaking result 
adding a new perspective to the study of non-classicality 
of the correlations in multipartite systems. By defin- 
ing entanglement operationally as those correlations that 
cannot be generated through local quantum operations 



assisted by classical communication, Werner noticed the 
existence of entangled states that do not violate any 
Bell's inequality. Surprisingly, Popescu showed that some 
non-separable-local states can be used to do quantum 
teleportation with fidelity greater than what is possible 
by classical means Q. 

Recall that, in quantum mechanics, the more general 
quantum operation describing changes between quantum 
states can be written, in the operator-sum representation, 
as Hi 



where Kj are the Kraus' operators acting on the state 
space of the system, satisfying the relation 



< 



with I being the identity operator in T-L. The equality in 
this relation is a sufficient condition for the map £ : H ^ 
"H to be a trace-preserving one, that is, Tr [£{p)] — 1. 

It turns out that the most general bipartite state that 
can be created via local quantum operations {£a and £b) 
and classical communication is of the form 



Sep 

Pab = 



(1) 
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where {pj} is a probability distribution (i.e. pj > 
J2jPj — 1) ^^'^ pf^^'' ^ 'Ha{b) is the reduced density 
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operator of subsystem A (B). States that can be written 
in the form shown in Eq. Q are said to be separable. 
All the states that are not separable are defined as being 
entangled (or non-separable). 

The non-local and non-separable aspects of the corre- 
lations in a composite system, as discussed above, are 
distinguishing features of the quantum realm. Never- 
theless, we believe that, in a more general fashion, it is 
the possibility of locally redistribute the correlations in 
a multipartite state that broadly defines its nature. In 
2008, Piani and co-authors proved the so called no-local- 
broadcasting theorem f9], showing that the correlations 
in a bipartite state pab can be redistributed via local 
quantum operations if and only if this state can be cast 
in the following form: 

PAi3=E^'^-fc|V'/)(V'/l®l<^f)(<^fL (2) 

where {pjk} is a probability distribution and 
iU't'k)}) orthonormal basis for subsystem A (B). 

By comparing this equation with Eq. ([T]) we see that 
p'^g is a subset of the separable states (which forms a 
convex set). Therefore, even separable states can present 
quantum correlations, in the sense that such correlations 
cannot be locally broadcast. This kind of non-classical 
correlation can be indicated, and in some sense quanti- 
fied, by a quantity called quantum discord (QD) p^[l6l|. 
that will be discussed in more detail in the next subsec- 
tion. 



B. Quantum Discord 

In this article we will focus on QD. This choice is 
motivated by the fact that this measure was the first 
one proposed in literature and was studied in most of 
the works presented so far. A detailed discussion about 
other measures of quantum correlations, as well as sev- 
eral other aspects of such measures, is presented in the 
reviews p^[l6l|. 

Before discussing QD it is instructive to briefiy review 
some aspects of classical information theory, in order to 
set the notation and define the relevant quantities. 

1. Some Concepts from Classieal Information Theory 

In classical information theory, the uncertainty about a 
random variable A is quantified by its Shannon's entropy 

m 

H{A) -.^ -Y^Palog^Pa, (3) 

a 

with Pa Pr(A = a) standing for the probability that 
A assuming the value a. The binary logarithm is used 
here such that information is measured in bits. 



Two random variables are said to be correlated if they 
share some information. Thus, by knowing the value of 
one of them, say B, we can acquire information about 
the other, A. The difference in the uncertainty about A 
before and after we know B, 

J{A:B):^H{A)~H{A\B), (4) 

is a measure of the correlations between the two random 
variables {A and B). This quantity is called mutual in- 
formation. In the above equation the conditional entropy 
reads 

H{A\B) := - ^Pa,b^Og2Pa\b, 
a, 6 

where Pa,b '■= Pr{A = a, B = b) and, from the definition 
for the conditional probability for A to assume the value 
a when B is equal to 6, Pa\b ■= Pa,b/Pb, it follows that 

H{A\B) = H{A,B) - H{B), 

with H{A,B) = ~J2a,bPa,b^og2Pa,b- Heucc an equiv- 
alent expression for the mutual information between A 
and B can be written as 

I{A:B) = H{A) + H{B) - H{A,B). (5) 

In classical information theory, the relation 

J{A:B) = I{A:B) 

is always valid. However, the extensions to the quan- 
tum realm of these two, classically equivalent, expres- 
sions for the mutual information differ in general. Such 
non-equivalence is at the heart of the quantum discord, 
that turns out to be a measure of quantum aspects of 
correlation, as will be discussed next. 



2. The Original Definition of Quantum Discord 

The uncertainty about a system S, described by a den- 
sity operator ps, is quantified in quantum information 
theory through its von Neumann's entropy (TT| : 

S{ps):=-TT{ps\og2Ps). (6) 

As a result, a direct extension for the mutual information 
shown in Eq. ([SJ can be obtained as 

T[pab) S{pa) + S{pb) - S{pab)- (7) 

This quantity is named quantum mutual information and 
it is regarded as a measure of the total correlations be- 
tween subsystems A and B when the state of the joint 
system is pab 

With the purpose of obtaining a quantum extension 
for the mutual information in Eq. Q , let us consider the 
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measurement of an observable represented by the follow- 
ing Hermitian operator 

OB = ^o,|nf)(nf|, 

3 

that is defined on the state space of subsystem B, Hb- 
If the system is initially in state pab, the value Oj is 
obtained with probability 

= Tt{Ia ® nfpAB) 

and the state of subsystem A, immediately after the mea- 
surement, reads 

pf = —TTBiU ^ nfpABU ® nf ). 

Pj 

In the above equations, {ll^} is a complete set of von 
Neumann's measurements on subsystem B satisfying the 
relations Y.j nf = Is and nf nf = SjkUf . 

The average quantum conditional entropy of subsys- 
tem A, given the measurement of the observable Ob on 
subsystem B, is given by 

SipA\B)^J2Pl^(pf)- 

3 

Thus, a quantum extension for J{A:B) can be defined as 

J{pab) ■■= S{pa) - S{pa\b)- (8) 

Equations (O and (|5]) for the quantum mutual infor- 
mation are not, in general, equivalent. The difference 
between them, 

VBipAB) ■=^iPAB) - ma.xJ{pAB), (9) 
Ob 

was named quantum discord [l3 | and is a measure of 
quantum correlation in bipartite systems. The maximum 
in Eq. is obtained through the measurements (on sub- 
system B) that provide the maximal information about 
A. It is worthwhile mentioning that another version for 
QD, Va, can be obtained if we measure an observable 
Oa on subsystem A. QD is an asymmetric quantity, in 
the sense that, in general, Va 7^ and has several 
interesting phy sical and information-theoretical interpre- 
tations [1^, [1^ . 

We also observe that if we use the mutual information 
as a measure of total correlations, one can verify that 
QD may be written as the difference between the mutual 
information of subsystem AB before and after a complete 
map of von Neumann's measurements is applied to one 
of the subsystems: 

T^BiPAB) = I{pab) - maxl(n^(p^s)), 

where Ii^{pAB) = Tji^A <E) nf )p^b(Ia <E) nf ). In this 
alternative definition, QD can be interpreted as a mea- 
sure of those correlations that are inevitably destroyed in 
the measurement process. 



C. NMR Systems in Quantum Information Science 

The actual relation between non-classical correlations 
and the advantages obtained in quantum information sci- 
ence (compared with classical protocols) is an important 
issue that must be thoroughly analysed in seeking for a 
better understanding about the source(s) of the quan- 
tum speedup. Nuclear Magnetic Resonance (NMR) has 
been one of the leading experimental platforms for the 
implementation of protocols and algorithms in quantum 
information science (QIS) (17i] and also for the simulation 
of quantum systems [1^. In NMR implementations, the 
qubits are encoded in nuclear spins that are manipulated 
through carefully designed sequences of radio-frequency 
pulses. The information about the system state is ob- 
tained directly from the transverse magnetizations, that 
are the natural observables in NMR experiments [s^]. 
The NMR density operator is typically highly mixed. In 
fact, Braunstein and co-authors showed that no entan- 
glement was generated in most NMR implementations of 
QIS protocols T^. Moreover, Vidal proved that a large 
amount of entanglement must be generated if a pure- 
state quantum computation is to provide an exponential 
speedup in information processing [20l |. The conjunction 
of these results led to a questioning about the possible 
classical nature of the NMR system and its QIS imple- 
mentations. On the other hand, the generation of en- 
tanglement was shown to be a necessary but not a suffi- 
cient condition for the existence of a gain, even in pure- 
state quantum computation [2l|. Moreover, the source 
of quantum speedup in the mixed-state quantum compu- 
tation scenario is even more subtle '22M25'|. In this case, 
QD seems to be the figure of merit for quantum advan- 
tage [2^, [l^l , although a definitive proof is still lacking. 
In this review we have a more modest goal. We deal only 
with the theoretical and experimental quantification and 
identification of quantum correlations in NMR systems 
at room temperature. We show the existence of quan- 
tumness in the correlations of highly mixed NMR states 
and study its peculiar dynamics under decoherence. 



II. QUANTIFICATION AND IDENTIFICATION 
OF QUANTUM CORRELATIONS 

A. Symmetric Quantum Discord 

As mentioned in the previous section, the original def- 
inition for QD is asymmetric with relation to the subsys- 
tem we choose to measure. To illustrate this fact, let us 
consider the following bipartite state 

piB-Y.p^\^f)(^^\^p^^ (10) 

i 

where {p,} is a probability distribution, {IV'i^)} G Ha is 
an orthonormal basis for subsystem A, and pf G Hb the 
reduced density operator of subsystem B. In this case. 
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the QD obtained by measuring subsystem A, 'Da{p'ab)^ 
vanishes while its alternative definition, 'Db{Pab)^ ^^^^ 
if and only if all pf commute, that is, if [pf^pf] = 
for all i and j. Moreover, numerical calculations indicate 
that whenever there exists an asymmetry in the state pab 
regarding the exchange between the subsystems or, more 
specifically, whenever S{pa) 7^ S{pb) it will result that 
T^a{pab) 7^ T^b{pab) • This asymmetry of quantum 
discord may seem weird in principle because one usually 
think about correlations as the information shared be- 
tween subsystems. We note that such interpretation still 
holds if one recognizes that the degree of quantumness 
of the subsystems shall limit the amount of information 
that the observers can extract via local measurements. 

The class of states in Eq. ([2]) is the only one for which 
T^a{Pab) — '^b{p'ab) = 0. a quantity that indicates 
and, in some sense, quantifies the quantumness of cor- 
relations in states that cannot be cast as p'^^g can be 
defined as (28i] 

V{pab) -^Apab) - maxI(n^^(pAB)), (11) 

where the complete map of local von Neumann's mea- 
surements reads 

n^^(p^s) := E(n/ ® ^k)PAB{i^f ® nf ), 

with nj = Its and n|n^ = (Jj^ni for A,B. This 
quantum correlation quantifier can be regarded as a sym- 
metric version of QD. Once state H'^^Ipab) is a classical 
one, we define the classical counterpart of I?, 

C{pab) ■■= maxZ(n^^(pAB)), (12) 
as a measure of the classical correlations in pab ■ 



Symmetric Quantum Discord for the Deviation 
Matrix 



In this subsection we show how to obtain the classical 
correlation (1121) and the symmetric quantum discord (I lip 
for systems described by density matrices of the form [33| 



By means of a Taylor expansion of ln(l/4-|-eAj) we obtain 

S{pab) = -7^Tr(pABlnpyis) 
in 2 

Iab , e(21n2-l) ^ 

In 2 ^'^^ 



2e2 
'ln2 



2-^Tr(Ap^,) 



(14) 



Once the reduced density matrices have the form 

with s = A, B, the same procedure can be used to com- 
pute 

S{Ps) = 1 - i^Tr(Ap2) + . . . . (15) 

Using (IT4I) and dTSl) it is possible to write the quantum 
mutual information in terms of the deviation matrix, up 
to second order in e, as 

I{pab) « ^[2Tr(Apis)-Tr(Api)-Tr(Ap|)]. (16) 

The measured state obtained from pab through a com- 
plete set of local-projective measurements reads 



U^''ipAB)^^-^ + eArjAB, 



with the measured deviation matrix given by 
AijAB :-n^^(ApAB) 

^^(n>nf)Ap^B(n/®nf) 

and Tr [U^^{Apab)] = 0. Thus, the same procedure 
utilized when obtaining leads to the following ex- 
pression for the mutual information of the measured state 

liU^^ipAB)) « ^[2Tr(A,7iB)-Tr(A,7i)-Tr(Ar;|)], 

(17) 

where Ar]A(B) = '^''^'B(A){AriAB)- The symmetric QD and 
the classical correlation are then computed as shown in 
Eqs. (fTT|) and HTH . respectively. 



PAB = + £ApAB, 



(13) 



C. Symmetric Quantum Discord for Two-Qubit 
States 



with e ^ 1 and ApAB being the traceless deviation ma- 
trix, Tr(A/9AB) = (Eq. is the typical density op- 
erator describing the state of NMR systems). In order 
to calculate the quantum mutual information we need to 
compute the von Neumann's entropy. For that purpose, 
let us use the eigen-decomposition of pab to write 

ln(p^B) = ^ln(l/4 + eA,)|A,)(A,|. 



Any two-qubit state can be brought, through local uni- 
tary transformations, to the following form 



1 



PAB 



Iab + ^{ajcrf (^iB + b^lA® af) 



,(18) 
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where {(t|} are the Pauh operators acting on Hs and 
Qj, hj, Cj S M. The measured state, modulo local unitary 
transformations, can be written as 



with 



Cjzfzf. 



The parameters 



with 



zf = 2 sin 6s cos 9s cos ( 
Z2 = 2 sin 9s cos 6*3 sin (j 
= 2cos2 9s - 1, 



G [0,7r/2] and (/.^ € [0,27r], 




FIG. 1. Schematic representation of the two systems A and 
B, initially prepared in a state pab, evolving under the action 
of local independent phase channels. 



D. Peculiar Dynamics of Correlations under 
Decoherence 



determine the measurement direction for the subsystem 
s = A, B. It seems that the maximization problem for 
obtaining the classical correlation cannot be solved in 
general. Nevertheless, for Bell-diagonal states. 



hd 
Pab 



(19) 



the correlations can be obtained analytically. In this case 



we have that Sip'f) = S^p'j^) = 1 and thus Eq.^ be- 



come 



with 



l + (-lp7, l + (-l)^7 

l0g2 ■ 



We note that the minimum of the measured state entropy 
is obtained by maximizing |7|. It results that 

I7I < k:= max(|ci|, |c2|, |c3|). 

Thus, an analytical expression for the symmetric QD for 
Bell-Diagonal states is obtained as [28|: 

1 

-Dip'^B) = E ^5" log2 4A5^ - C{p%), (20) 

j,k=a 



with 



1 1 

C{p'^b) = ^ E[1 + (-1)'^] l°g2[l + i~^y^] (21) 



and X'^i = [l + (-l)^ci - {-iy+''c2 + (-l)Jc3]/4. 



We shall use the analytical formulas for the symmetric 
QD (PO)) and for its classical counterpart (PT|) in Bell- 
diagonal states to study the dynamics of these correla- 
tions under the influence of local independent Markovian 
environments that inject phase noise into the system, as 
schematically depicted in Fig. [T] 29] . Under such con- 
ditions, the two subsystems A and i?, that are initially 
prepared in a Bell-diagonal state ([T9)) . evolve via the fol- 
lowing quantum operation [sO] 

PABiP) = SpDip'^B) 

= E(^/ ® Kk)PABiKf (22) 

"1+7 a-p- 

_1 l-ja + /3 

^4 a + 13 1-j 

a~/3 1 + 7. 



(23) 



All matrices of the article are represented in the stan- 
dard computational basis {|00), |01), |10), |11)}, where 
:— \i) (8) |j) and {|0), |1)} are the eigenstates of the 
Pauli matrix CTz. The Kraus' operators appearing in the 
operator-sum representation (|22p for the phase-damping 
channel are 



A'" — 



1 - 

Ps 

2 



1 
1 



1 

-1 



(24) 



with s — A, B and pa = Pb '■= P being the parametrized 
times. This means that we considered identical environ- 
ments for both subsystems. Besides, we have defined 

(l-p)'ci, 

(1-P)'C2, 
C3- 



0.00 

0.0 0.2 0.4 0.6 0.8 1.0 
P 

FIG. 2. Total correlation (blue dotted line), classical correla- 
tion (red dashed line), and quantum discord (gray continuous 
line) for a Bell-diagonal initial state evolving under local in- 
dependent phase-damping channels. For this figure, the pa- 
rameters in the initial state p9p were chosen as: ci = 0.06, 
C2 = 0.3, and C3 — 0.33. For this state the classical corre- 
lation is not affected by the environment while the quantum 
correlation decays monotonically. 

We see that the evolved density matrix ([23)) is in Bell- 
diagonal form p^ . Thus, the symmetric QD V [p'^^^ip)] 
and the classical correlation C [p'^sip)] given by ([20)1 
and (|21l) . respectively, with 

K = max(|a|, |^|, I7I) 

and 

= [1 + i-iya i-iy+'^p + i-iy^]/4. 



From the analytical results we readily identify three 
general and different kinds of dynamic behaviour of the 
correlations under decoherence: 

: (i) If jcal > |ci|, |c2| we have k = jcsl and the classical 
correlation does not depend on the parametrized 
time p. As the quantum correlation decays mono- 
tonically going to zero in the asymptotic state limit, 
the classical correlation is equal to the mutual in- 
formation of this state, i.e. 

c [p'Mp)] = c [p'Mp = o)]=i [p'Mp = 1)] ■ 

This kind of dynamical behaviour of the correla- 
tions is exemplified in Fig. [5] for the initial state 
p^iB with ci = 0.06, C2 = 0.3, and C3 = 0.33. 



0.0 0.2 0.4 0.6 0.8 1.0 
P 

FIG. 3. Total correlation (blue dotted line), classical correla- 
tion (red dashed line), and quantum discord (gray continuous 
line) for the a Bell-diagonal initial state evolving under lo- 
cal independent phase-damping channels. For this figure, the 
parameters in the initial state H19p were chosen as: ci = 1, 
C2 — —0.6, and C3 = 0.6. For this state the sudden change 
occurs at psc ~ 0.22 and the quantum discord (classical corre- 
lation) remains constant (decays monotonically) for p < 0.22 
with the opposite scenario taking place for p > 0.22. 



Classical correlation decays monotonically until a 
specific parametrized time 

P.c = i-\ — rrrn^^ (25) 

y max(|ci|, \C2\) 

and remains constant form this time afterwards. 
For p < Psc we have k = (1 —p)^ max(|ci|, |c2|) and 
C decays monotonically. On the other hand, for 
P ^ Psc we have k = jcaj and C is constant: 

c[p%{p>psc)] ^lip'Up-^)], 

and the decay rate of V changes abruptly in p = 
Psc- The dynamical behaviour of correlations with 
a sudden change in its decay rate is exemplified 
in Fig. [3] for the initial state p'^g with ci — 1, 
C2 — —0.6 and C3 = 0.6. For this example, not only 
the classical correlation remains constant in a cer- 
tain time interval but also the quantum correlation 
can be immune against decoherence. The regimes 
where V = constant and C = constant were named 
classical and quantum decoherence regimes in Ref. 




: (iii) Finally, if C3 = both correlations C and V decay 

: (ii) If |ci| > |c2|, Icsl or |c2| > |ci|, jcal and jcal ^ we monotonically. We exemplify this kind of dynamic 

have a peculiar dynamics of the correlations under behaviour in Fig. Hfor a Bell-diagonal initial state 

decoherence with a sudden change in behaviour. with ci = C2 = 0.25 and C3 = 0.0. 
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FIG. 4. Total correlation (blue dotted line), classical corre- 
lation (red dashed line), and quantum discord (gray contin- 
uous line) for the a Bell-diagonal initial state evolving un- 
der local independent phase-damping channels. For this fig- 
ure, the parameters in the initial state (|19[l were chosen as: 
Ci = C2 = 0.25 and C3 — 0.0. For this state all correlations 
decay monotonously with time. 



E. Classicality Witness 

The set composed by separable states of the form ([T|) 
is convex and hence, there exists a linear hermitian op- 
erator, called entanglement witness, that can be used to 
distinguish separable from entangled states. It is natural 
to ask whether the same approach can be utilized to iden- 
tify classical states. It turns out that the set formed by 
states with null QD is not convex and, as a consequence, 
linear witnesses cannot be applied any more. To verify 
this assertion, let us consider a linear hermitian operator 
Wi whose mean value is non-negative for classical states, 
i.e. 



Pab 



Tr(iy,p5fB) > 0, 



and (VFi)p^B < for quantum-correlated states. Thus 
we note that for the separable states (H]), 



Therefore linear witnesses cannot be applied to identify 
separable discordant states. Next we introduce a non- 
linear classicality witness providing a sufficient condition 
for absence of quantumness in the correlations of two- 
qubit states of the form (IT51) [32*1 . First we consider the 
following set of hermitian operators 



04 = z-aA 



where j = 1,2,3 and z,w eM.'^ {\z\ — \w\ = 1) should be 
randomly chosen. Let us consider a relation among the 
averages of these operators as follows 



3 4 

E E 

i=l i=l+l 



(26) 



where |a;| is the absolute value of x. We see that Wp^^ = 
if and only if the mean value of at least three of the 
four observables defined above are zero. As 

{Oi)pAB = for i = 1,2,3, 
{04:)pAB = z-d+w -b, 

with a ~ (oi, 02, 03) and b — (61, &2, ^3), the only way to 
warrant that Wpj^^ = is if pab is of the type 

Xi — -ri^AB + Ciaf (g) erf) for i = 1 or i = 2 or i = 3, 



X4 



E(« 



'Ib +hjlA(E>af) 



It turns out that all theses states can be cast as those in 
Eq. ^ and, therefore, do not possess quantum correla- 



tions. As a consequence, it follows that Wp 



is a 



sufhcient condition for pab to be classically correlated or 
to have no correlations at all. 

For Bell-diagonal states, Wp^g = is also a necessary 
condition for the absence of quantum discord in the sys- 
tem. This result is obtained by noting that p^LB being 
classically correlated implies that it must have the form 



(I 



AB + Cta, 



)o-f )/4 with i 



1 or z = 2 or i = 3, implying 
that Wp^g = 0. 

For the experimental implementation of the classicality 
witness in the NMR context [? ], that will be pre- 
sented bellow, it is useful to rewrite the witness W, 



PAB 



in terms of the qubits magnetizations, that are the nat- 
ural observables accessed in NMR experiments. First we 
observe that 



{0^ 



{(7^ ® Is) 



ith 



where 



= UA^B[RnM)PABRUO,)]UA 



^B, 



Ua^b = |0)(0|®Is + |l)(l|®(Tf 

is the controUed-NOT gate with qubit A as control and 
i?„^(6»i) = R^^{Oi) (g) R^^iOi) with Rf,^ie,) being a local 
rotation on qubit s ^ A, B by a.n angle 9i {9i = 0, 6*2 = 
^3 = 7r/2) in direction hi (712 = y, = z). 



III. QUANTUM DISCORD IN NMR SYSTEMS 

A. The System 

In order to experimentally follow the dynamic of cor- 
relations under decoherence and to implement the clas- 
sicality witness discussed in the preceding section, it is 
necessary to prepare, manipulate, and measure the state 
of a two-qubit system in the laboratory. In the NMR sce- 
nario, two-qubit systems can be achieved using samples 
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presenting either two J coupled nuclear spins- 1/2 (where 
J means scalar spin-spin coupling) or one quadrupolar 
spin-3/2 system in the presence of a local electric field 
gradient (named as quadrupolar spin systems). There 
are many examples of two dipolar J coupled spins- 1/2 
in NMR, for instance, and ^^C nuclei in chloro- 
form (CHCI3) or and "^^P nuclei in phosphoric acid 
(H3PO4) [ITi]. Two-qubit quadrupolar systems used so 
far in QIP are usually comprised by spin-3/2 nuclei in 
single crystals [s^ or in lyotropic liquid crystals [13, ■ 
^^Na and ^Li nuclei in lyotropic liquid crystals based on 
Sodium Dodecyl Sulfate (SDS) [37] and Lithium Tetraflu- 
oroborate (LiBF4) [39|], as well as ^'^Na in Sodium Nitrate 
(NaNOa) single crystals [s^ are some examples of these 
quadrupolar systems. 

All spin systems cited above are equally good represen- 
tations of two-qubit sytems. However, the nuclear spin 
interactions that drive their quantum evolution are dis- 
tinct, leading to significant differences in the techniques 
used for state preparation, manipulation and read-out. 
Besides, as a result of the specific characteristics of each 
spin interaction, there are unique features in their deco- 
herence behaviour. Hence, in the next sections we shall 
discuss each system separately. 




Frequency Frequency 



FIG. 5. Top: Schematic representation of Chloroform 
molecule (HCCI3). The red (blue) circles depict perturbation 
and response of Carbon (Hydrogen) nucleus at its respective 
resonance frequency. Bottom: The normalized equilibrium 
spectra of Carbon (Red) and Hydrogen (Blue) nuclei. 



1. Two Dipolar J Coupled Nuclear Spin-1/2 Systems - 
Liquid Sample 

A schematic representation of a ^^C enriched chloro- 
form molecule (CHCI3) is depicted in Fig. [SJ As already 
mentioned, the two-qubits are comprised by the nuclear 
^H — ^'^ C spin pair in the molecule. However, due to sen- 
sitivity reasons, one need to have many replicas of these 
two-qubit units, i.e, the system need to be constituted 
by an ensemble of spin pairs where each spin in the pair 
interacts only with its counterpart. This is achieved by 
diluting the ^^C enriched CHCI3 molecules in a deuter- 
ated solvent. Since the deuteration and the low natural 
abundance of ^"^C make the solvent molecules magneti- 
cally inert, in the NMR point of view this system can be 
well described as an ensemble of isolated spin pairs, i.e, 
many replicas of two-qubit units. Thus, all the results 
reported here were obtained using a sample prepared by 
dissolving 100 mg of 99% i^C-labeled CHCI3 in 0.2 ml of 
99.8% Acetone-d6, and placing this in a Wildmad Lab- 
Glass 5 mm tube. Both samples were provided by the 
Cambridge Isotope Laboratories - Inc. NMR experiments 
were performed at 25° C using a Varian 500 MHz Spec- 
trometer and a 5-mm double-resonance probe tuned to 
IH and 13C nuclei. 

In NMR experiments the nuclear spins are placed in a 
strong static magnetic field Bq, whose direction is taken 
to be the positive z. Spins manipulation is achieved by 
applying time dependent radio-frequency (rf) fields on 
both spins, which interact to each other and with their 
local environments. Thus, in a double rotating frame 
with rf frequencies and oj^f, the Hamiltonian of a 
two coupled spins is represented by 

n ^-{ujH - Olf - (c^c - c.f/)lf + 2^Jlf 
(If cos(^^ + I«sin^«) 
+uj^ (I^ cos + sin (^^) + Ueuv (t), (27) 

where (I?) the spin angular momentum operator 
in the u,v — x,y, z direction for ^H (^'^C); (p^ 
LUi (wp) define phase and power, respectively, of the 
radio-frequency field for the ^H (^'^C) nuclei. The first 
two terms describe the free precession of ^H and ^■^C nu- 
clei about Bq, with Larmor frequencies 

— « 500 MHz and — w 125 MHz 
27r 27r 

The third term describes the scalar spin-spin coupling 
(also refereed as J coupling in the NMR literature) of 
frequency 

J « 215.1 Hz. 

The fourth and fifth terms are the external rf fields ap- 
plied to manipulate the ^H and ^^C nuclear spins, respec- 
tively. HEnvit) represents time dependent fields result- 
ing from the random fluctuating in the interactions be- 
tween the spins and their environment. This term leads 



to spin relaxation and decoherence and includes inter- 
actions with the chlorine nuclei as well as higher order 
terms in the spin-spin coupling. 

The density operator for any quantum system in con- 
tact with a thermal bath can be represented as 



P 



Trexp(-/3-H)' 



(28) 



where /3 = l/ksT with ks being Boltzmann's constant 
and T the absolute temperature. For two-qubit NMR 
systems at room temperature, the thermal energy is 
much greater than the magnetic energy: 



^ksT 
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where 7 is the gyromagnetic ratio. In this case, it is 
possible to write the density operator in the so called 
high-temperature expansion: 



lab 
T 



eAp, 



(29) 



with Ap being the deviation matrix. We observe that 
Ap is the part of the system density matrix which is 
actually manipulated (via rf pulses) and accessed in NMR 
experiments. 

Two unitary transformations are frequently used in 
NMR experiments: radio-frequency pulses and free evo- 
lutions of the spin systems under the coupling and static 
fields. The rf pulses are characterized by some intrin- 
sic parameters as amplitude (power) , duration, frequency 
and phase. By proper set of these parameters, it is pos- 
sible to induce nuclear spin rotations about any axis and 
with any rotation angle. For example, a non-selective 
spin rotation (also known in the NMR jargon as hard 
pulses) of 7r/2 can be performed by a rf pulse with du- 
ration of 7.4 /is. Using these time values we acquire the 
equilibrium spectrum (NMR spectrum acquired after ap- 
plying the rf pulses to the thermal equilibrium state) for 
each nuclear species (see Fig. [5]). Analogously, the free 
evolution corresponds to precession of the nuclear spins 
around Bq. 



2. Nuclear Quadrupolar System - Liquid Cristal Sample 

As already mentioned, ^^Na (/ = 3/2) in a SDS 
molecule of a lyotropic liquid crystal is a good repre- 
sentation of a two-qubit quadrupolar spin system. As 
depicted in Fig. [51 this molecule has a chemical struc- 
ture were the ^^^Na nucleus is far from other NMR ac- 
tive high abundant nuclei, mostly ^H, so the magnetic 
dipolar interaction of ^^^Na with other nuclei in the SDS 
molecule can be neglected. Besides, since the liquid crys- 
tal is prepared with heavy water (D2O), the dipolar in- 
teraction with the nuclear spins of water is also week 
(deuterium is a low 7 nucleus). Thus, in the SDS D2O 
based liquid crystal, the interaction between the ^'^Na 




FIG. 6. Left-top: Labels of atomic compounds of the 
molecules. Left-bottom: The normalized equilibrium spec- 
trum of the Sodium nucleus. Right: Cartoon of the Sodium 
Dodecil Sulfate molecule (SDS). The green ellipses depict the 
perturbation and response of Sodium nucleus at its respective 
resonance frequencies. 



quadrupole moment and the electric field gradient pro- 
duced by the electrical charges in the vicinity (usually 
named as quadrupolar coupling) can be considered as 
the only internal spin interaction that affects the quan- 
tum system evolution. However, the local electric field 
gradient is defined by the charges configuration in the 
^■^Na nucleus proximity and to produce a system were 
each ^'^Na environment is replicated along the sample, all 
SDS molecules needs to have similar orientations. Fortu- 
nately, the strong static magnetic field applied in NMR 
experiments induces the alignment of the SDS molecules 
in the lyotropic liquid crystal, naturally producing a sit- 
uation were the local electric field gradients are identical 
at each ^'^Na site. Furthermore, the anysotropic inter- 
nal motions of the SDS molecules reduces the strength of 
the quadrupolar coupling to smaller values, facilitating 
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the manipulation of the spin system by rf pulses. Under 
such conditions, in a rotating frame with frequency oJrf, 
the ^^Na nuclear spin Hamiltonian can be described by: 

H = - K - iOrf) I. + ^ (31' - I') 

+UJ1 {I^ cos ^ + Iy sin <f ) + n%^^ (t) , (30) 

where loq is the strength of the quadrupolar coupling in 
frequency units and cjl is the Larmor frequency {\ujl\ ^ 
IwqI). The spin angular momentum operators are repre- 
sented by its M-component I„ {u — x, y, z) and its square 
modulus I^. The first term of Eq. ([50]) describes the 
Zeeman interaction between a static magnetic field and 
the nuclear spin, while the second term accounts for the 
static first order quadrupolar interaction. As in Eq. (|27)) , 
the third and fourth terms correspond to an externally 
applied rf field and the time dependent coupling of the 
^^Na nuclear spin with its environment, respectively. In 
this case, Ti-^nvi^) results mainly from the motion in- 
duced fluctuations in the local electric field gradients with 
smaller contributions from the weak dipolar interactions 
between the ^^Na and other NMR active spins in the 
surroundings [iol. l42j|. 

The experimental results presented here were obtained 
using a liquid crystal sample prepared with 20.9% of SDS 
(95% of purity), 3.7% of decanol, and 75.4% deuterium 
oxide, following the procedure described in Ref. [43]. The 
^■^Na NMR experiments were performed in a 400 MHz - 
VARIAN INOVA spectrometer using a 7 mm solid-state 
NMR probe head at 26°C. In the left-bottom of Fig. Elwe 
show the ^"^Na equilibrium spectrum in our SDS liquid 
crystal sample. The quadrupole coupling frequency is 
directly obtained from the separation between the central 
and the satellite lines [4^, |45|], which in this particular 
case gives 



"^^Ib kHz. 
27r 



Once the general features of the two main two-qubit 
NMR systems were presented, we now wish to discuss 
how these systems can be used for QIP. To do so, one 
needs to show how to perform the three main QIP steps, 
i.e state preparation, manipulation and read-out. This 
will be the focus of the next three subsections. 



B. State Preparation 

Most of the NMR QIP implementations relies on effec- 
tive pure states, usually refereed as pseudo-pure states. 
Among others |46l - t48| . the main procedures for pseudo- 
pure state preparation are the so called spatial [47] and 
temporal averaging [46j ]. which we shall discuss in this 
section. 

The general idea behind the spatial averaging is, start- 
ing from the thermal equilibrium, to apply a set of rf 
and magnetic field gradient pulses to manipulate the spin 



populations and coherences (diagonal and off-diagonal el- 
ements of the density matrix in the computational basis, 
respectively) in order to produce the desired quantum 
state. While the rf pulse acts changing the density ma- 
trix populations and coherences through specific spin ro- 
tations, the field gradients are used to clean-up undesired 
coherences. Indeed, the application of the gradient pulses 
characterizes the spacial averaging, since they produce a 
magnetic field distribution along the sample that induces 
a random phase distribution in the individual quantum 
coherences, so they do not contribute to the final average 
state. For instance, the quantum state ]11) of the compu- 
tational basis can be prepared by applying the following 
pulse sequence to the thermal equilibrium state: 



' ~V 



' V 

k2j) 



ii.JJ 



(f) 



H,C 



where {d)^''~^ is a rotation by an angle 6 in d direction ap- 
plied simultaneously on ^H and ^'^C nuclear spins. (1/ nJ) 
represents the duration of a free evolution. Last {Gz} 
states for a magnetic field gradient in z direction. Note 
that, in this case, the specific rotations are implemented 
by hard rf pulses with proper phases (sub indices) and 
amplitudes (rotation angle between parenthesis) applied 
to ^H and -^^C (super indices) nuclear spins. 

The pulse sequence (rf pulses, free evolutions, and 
magnetic field gradients) for preparation of the state ]11) 
was numerically simulated and the real and imaginary 
parts of the deviation matrix are shown in the bar plots 
in the left part of Fig. [71 The experimental results, and 
corresponding pulse sequence, are displayed on the right 
part of Fig. [T] To obtain the experimental deviation ma- 
trix, the quantum state tomography procedure, which 
will be described in Sec. IIII C[ was used. 

State preparation by spatial averaging can also be 
achieved replacing the hard rf pulses and free evolutions 
by numerically optimized low power pulses (known in 
the NMR terminology as soft pulses) followed by mag- 
netic gradient pulses. Due to the excellent control of the 
NMR parameters that dictate the system evolution, the 
NMR Hamiltonian can be parametrized in terms of the 
rf pulses parameters (power wi, duration t and phase ip). 
This allow using numerical optimization routines to seek 
for the set of pulses, and the corresponding parameters, 
that brings the system from the thermal equilibrium to a 
state whose the diagonal element of the deviation matrix 
are the same as those of the target state. Thus, the sub- 
sequent application of the field gradient pulse vanishes all 
off-diagonal deviation matrix elements, producing the de- 
sired state. Two numerical procedures named as GRAPE 
(GRAdient Ascent Pulse Engineering) [H [H^ and SMP 
(Strongly Modulated Pulses) 51] are commonly used for 
pulse optimization. An example of a SMP pulse sequence 
used to implement an specific initial quantum state is pre- 
sented in Table H] The first column indicates the order of 
application for a certain set of parameters. The second 
(fifth) and third (sixth) columns show the numerical val- 
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Numerical procedure 

Real 



Experimental procedure 



Real 





Imaginary 



Imaginary 





trix shown in the left part of Fig. |51 An experiment 
using the parameters values shown in Table |T] was also 
performed. The resulting deviation matrix, also obtained 
using quantum state tomography, is shown in the right 
part of Fig. [S] 



Numerical procedure 

Real 




Experimental procedure 

Real 



^^^^^ 




Imaginary 



Imaginary 



FIG. 7. Bar representation of the prepared deviation ma- 
trix. These states were produced by using hard pulses, free 
evolutions, and magnetic field gradients. The top (bottom) 
part of the figure corresponds to the real (imaginary) part of 
the deviation matrix while the left (right) part of the figure 
corresponds to the simulated (experimental) results for the 
prepared initial quantum state. 



SMP 


Parameters 


Pulse 


Hydrogen 


time 


Carbon 


Step 


^? 




(ms) 






1 


1092.6 


5.94 


0.518 


1962.9 


1.88 


2 


1679.8 


5.73 


1.546 


1276.7 


5.26 


3 


112.0 


1.64 


0.226 


619.3 


2.55 


4 


1998.0 


2.63 


1.130 


868.9 


3.91 


5 


211.3 


4.84 


1.996 


835.8 


5.46 



TABLE I. Numerical values of the power, duration and phase 
of the rf pulses utilized to implement the SMP technique for 
quantum state preparation. These parameters were used to 
produce diagonal initial quantum states (see, e.g.. Fig. [TT]). 
The SMP technique is one of the steps used to produce suit- 
able quantum states for the study of quantum discord in NMR 
systems. For instance, using these values for the parameters 
characterizing the rf pulses we produced a Bell-diagonal state 
satisfying the following inequality: |ci|, \c2\ > |c3|. Similar 
procedures were followed to obtain other states of interest. 



ues of the power and phase of the rf pulses applied to the 
(^^C) nuclei. The duration of each rf pulse is shown 
in the fourth column. As already mentioned, after the 
RF pulse a magnetic field gradient is applied to wash out 
the off-diagonal coherences in the deviation matrix. 

The application of the pulse sequence showed in Table 
m followed by the gradient pulse, to the thermal equi- 
librium state was numerically simulated being the corre- 
sponding real and imaginary parts of the deviation ma- 





FIG. 8. Bar representation of the prepared deviation matrix. 
These states were produced by using the SMP technique. The 
top (bottom) part of the figure corresponds to the real (imagi- 
nary) part of the deviation matrix while the left (right) part of 
the figure corresponds to the simulated (experimental) results 
for the prepared initial quantum state. 

The main advantage of usiirg SMP or GRAPE pulses 
followed by magiretic field gradient pulses for state prepa- 
ration via spatial averaging is their versatility in terms of 
producing optimized pulses that leads to any desired spin 
populations, which makes the method more general than 
the use of hard pulses, where a specific pulse sequence 
has to be found for each desired state. 

The preparation of a quantum state via temporal aver- 
aging involves the realization of a specific number of time 
uncorrelated experimeirts. However, each experiment is 
designed in such a way that the resulting state obtained 
after summing up the result of all experiments, is the 
desired state. As in the case of spatial averaging, both 
the hard [i^ and soft pulses [H, [5l| approaches can be 
used. One advantage of this preparation scheme is the 
possibility of preparing non-diagonal states, that can be 
prepared directly from the thermal equilibrium, since no 
magnetic gradient pulses are applied. 

An example of a state prepared using the temporal 
averaging scheme is shown in [9] This state was imple- 
mented using ^^Na quadrupolar nuclei (spin 3/2) in a 
lyotropic liquid crystal sample as described in the exper- 
imcirtal sectioir. A SMP pulse was also used in the state 
preparation, but in this case the parairreter values were 
optimized based on the Hamiltonian of Eq. dSH]) . Dif- 
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Group 


SMP 


Parameters 




Pulse 


of Sodium nuclei 


Set 


Step 






time (/is) 




1 


37485.4 


6.208 


15.11 




2 


11624.8 


4.243 


36.53 


A 


3 


31867.5 


3.293 


4.15 




4 


38906.5 


6.027 


47.99 




5 


38681.0 


1.927 


31.62 




1 


12063.7 


3.919 


7.95 




2 


15399.8 


3.581 


33.74 


B 


3 


38721.2 


5.845 


30.47 




4 


24148.3 


5.117 


3.38 




5 


38916.9 


3.762 


48.08 




1 


39208.2 


1.788 


38.82 




2 


7589.9 


4.984 


27.38 


C 


3 


19817.7 


0.889 


9.27 




4 


38915.2 


1.682 


25.45 




5 


38755.6 


1.200 


21.98 




1 


13372.2 


4.333 


15.24 




2 


2192.1 


1.506 


37.89 


D 


3 


39094.7 


4.433 


29.35 




4 


11495.7 


3.294 


10.49 




5 


18867.5 


0.213 


36.92 



Numerical procedure 

Real 



TABLE II. Numerical values of power, duration, and phase 
of the rf pulses used to prepare quantum states by imple- 
menting the SMP technique and the time average procedure. 
These pulse sequences were used to produce the diagonal ini- 
tial quantum state (see Fig. Ilip . The SMP technique is one 
of the steps needed to produce appropriate quantum states to 
study quantum discord in NMR systems. Using these param- 
eter values we have produced quantum states that satisfy the 
following inequality: jci| > |c2|, |c3| 7^ 0. 



ferent experiments are labelled by the capital letters A, 
B, C, and D, each one having five sets of values for the 
rf pulses parameters (wi, t and v?), see Table HIl Each 
pulse was applied separately to the thermal equilibrium 
state, resulting in transformed deviation matrices, whose 
combination (sum) is the final Ap. The final simulated 
Ap is shown in the bar plot in the left of Fig. [S] The 
corresponding experimental Ap, obtained using same pa- 
rameter values of Table HIl is shown in the right side of 

Fig. m 



C. Quantum State Tomography 

In general, the Quantum state tomography (QST) uses 
the observables of the quantum system to quantitatively 
measure each element of the density matrix that corre- 
spond to the state of the system. In NMR systems these 
observables are comprised by spin magnetizations, which 
depend on a combination of the density matrix that 




Experimental procedure 

Real 




Imaginary 



Imaginary 





FIG. 9. Bar representation of the prepared deviation ma- 
trix. These states were produced by using hard pulses, free 
evolutions and magnetic field gradients. The top (bottom) 
part of the figure corresponds to the real (imaginary) part of 
the deviation matrix while the left (right) part of the figure 
corresponds to the simulated (experimental) results for the 
prepared initial quantum state. 



correspond to single quantum transitions (Am = ±1)- 
Thus, it is usually stated that only those elements lo- 
cated at certain reading position in the density matrix 
are experimentally accessed directly. Thus, in order to 
obtain the elements out of the reading positions it is nec- 
essary to apply a particular scheme of rotations to bring 
the desired elements to the reading positions in a con- 
trolled way. By doing so, the spin magnetizations be- 
come dependent on these "non directly detectable" el- 
ements, which can then be obtained after proper data 
processing. Thus, the art behind QST is to find a min- 
imum number of rotations that allow mapping the full 
sytem density matrix (actually the deviation matrix in 
the NMR context). Indeed, there are many QST pro- 
cedures in NMR QIP either using single or combined 
spin rotationsjlljIlHl, global rotations of the spin system 
(53! , transition selective excitations [s^, [4l| , among oth- 
ers. Generally speaking, each QST method is more ap- 
propriated for a specific spin system, for example, while 
the methods based on a set of single spin rotations are 
usually applied for spin 1/2 systems, the methods based 
on transition selective pulses or global spin rotations are 
more appropriated for quadrupolar nuclei. 

As a specific example, we will briefiy discuss here the 
QST technique introduced by Long et al. in Ref. [sj], 
a method widely used for QST in two-coupled spin 1/2 
system, which was also adapted for systems with higher 
number of spins [s^ [s^. In this method, the QST is 
performed in 9 steps and, before each step, the system 
needs to be prepared in the same state, whose the de- 



viation matrix need to be characterized. The first step 
consists in acquiring the NMR signal of both nuclei (free 
induction decay - FID) in the absence of radio-frequency 
pulses. After the Fourier transformation of the FID, the 
line intensities in the resulting spectrum, labelled as II 
in Fig. [TUJ are measured and recorded. In each of the fol- 
lowing steps, rf pulses with specific phases, as described 
in Table IIIIl are applied to each or both nuclei and the 
same reading procedure is executed. 

Fig. \TU\ shows in red (blue) the Carbon (Hydrogen) 
spectrum obtained after each tomography step, executed 
after preparing a Bell-diagonal state in the CHCI3 two 
qubit spin system. The spectra are labelled as II, XX, 
IX, lY, XI, YI, XY, YX, YY, with I representing no 
pulse and X (Y) a 7r/2 pulse applied in positive x (y) di- 
rection. In this case, since each spectrum has two spec- 
tral lines with real and imaginary components, 72 line 
intensities are recorded. Each of this line intensities are 
associated with one or more elements of the deviation 
matrix, so to obtain the full deviation matrix one needs 
to know the relations between the line intensities (NMR 
read-outs) and the deviation matrix elements. Indeed, 
the rf pulse phases on IIIIl are such that all the elements 
are accessible, with some redundancy for error minimiza- 
tion. 

To be more specific, let us consider the following devi- 
ation matrix, 

ai a2 + iaii 03 -I- iai2 04 -I- zaia 

03 — iai2 ag — «ai4 as ag -I- laie 

04 — zai3 a? ~ *ai5 0,9 ~ *ai6 0,10 

(31) 

which is described by sixteen unknown variables. The ro- 
tation imposed to the spin system by the rf pulses in each 
step n of the tomographic process can be represented by 
the operator Un(X,Y). Note that the analytical form 
of U„(X, Y) is known for a specified pulse duration and 



Pulse order 


Carbon 


Hydrogen 


1st 


I 


I 


2nd 


X 


X 


3rd 


I 


X 


4th 


I 


Y 


5th 


X 


I 


6th 


Y 


I 


7th 


X 


Y 


8th 


Y 


X 


9th 


Y 


Y 



TABLE III. An ordered pulse sequence applied simultane- 
ously in both nuclear species. X = (f)_,_^ represents a 
radio- frequency pulse of 7r/2 in positive x-direction. Y = 
represents a radio- frequency pulse of tv/2 in positive 
j/-direction. I represents absence of radio-frequency pulse. 
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FIG. 10. Set of spectra of a two-spin system encoded in the 
chloroform molecule generated in the implementation of QST 
pulse sequence detailed in Table Hill Left (right) spectra cor- 
responds to the Carbon (Hydrogen) nuclear spin. The bar 
plots in the middle corresponds to the real (top) and imag- 
inary (bottom) parts of the experimentally deviation matrix 
obtained through quantum state tomography. 

phase. Hence, after the n-th rotation, the transformed 
deviation matrix is represented by 

Ap„ = U„ (X, Y) ApUt (X, Y) . (32) 

Thus, after the nth transformation, the real (Mr) and 
imaginary (My) parts of the NMR magnetization for each 
nuclear spin A,B, become 

m4^„ = Tr [Ap„ (I,,, <E> I)] , (33) 
Mi^^^„ = Tr[Ap„(I®I,,,)], (34) 

From Eq. [MJ one can see that the magnetization val- 
ues are linearly related with the Ap elements, i.e., each 
n value in Eg 1341 represents a linear equation having the 
flj 's Ap elements as independent variables, the magneti- 
zations M^-^^ as dependent variables and the elements 
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of the operator Un(X,Y) as coefBcients. Note that the 
magnetization values are experimentaUy accessed and di- 
rectly related to the line intensities and the Un(X, Y) 
elements are known from the rf pulse parameters. Thus, 
Eqs. [331 can be rewritten as 



bk — Xhia 



kjU-j, 



(35) 



where bk with k = 1, 73 state for the fc-th recorded line 
intensity refer to the Ap elements before the transforma- 
tion and j = 1, 16. The elements Xkj are obtained 
from the pulse parameters applied in each QST step. 

The solution of equations (j35p for the aj variables en- 
able us to reconstruct the deviation matrix. An example 
of this reconstruction procedure is presented in the bar- 
graph of Fig. [TUl 

The QST methods for quadrupolar nuclei follows an 
analogous idea, i.e., obtaining a series of read-outs af- 
ter applying specific rf pulses, recording the correspond- 
ing line intensities and solving a linear equation system 
to obtain the deviation matrix elements. However, the 
specificities of the method are somewhat more intricate 
and will not be presented here. More information can be 
found in Ref. 1531. 



D. Experimental Measurments of Quantum 
Correlations in NMR systems 

In this section, experimental investigations regarding 
the quantification and identification of quantum correla- 
tions in NMR systems at room temperature will be pre- 
sented. At first, we verify that the peculiar dynamics of 
quantum discord under decoherence, theoretically pre- 
dicted for phase-noise channels [29j , can be present even 
under the influence of an additional thermal environment 
(33 [. Next, we present the first experimental implemen- 
tation of a classicality witness [32], which was used to 
identify the nature of the correlations in NMR systems 
without the necessity of a full QST f35| . 



1. Experimental Dynamics of Correlations 

Here we present an experimental verification of the 
sudden change phenomenon of quantum discord under 
decoherence, whose theory was discussed in Section III Dl 
Although, in principle, such kind of phenomenon may 
take place in the quadrupolar system in a liquid-crystal 
sample, introduced Section UlI A 21 33] . it is very difficult 
to observe it experimentally. The reason is that in this 
system the phase noise environment is global [5^, and 
the transversal relaxation time is very short, particularly 
when compared with that for the two nuclear spins-1/2 
in a liquid sample (see Section IlII All) . Thus, for experi- 
mentally demonstrating the peculiar behaviour of corre- 
lations under decoherence, a two-qubit system encoded 
in the ^H and "'^'^C nuclear spins in a CHCI3 molecules. 



as discussed in section lTlI A II was used. As already men- 
tioned, in this system the relaxation process that causes 
phase decoherence and energy dissipation is due, mainly, 
to internal molecular or atomic motions that leads to 
random fiuctuations in the electromagnetic field in which 
the qubits are immersed. The decoherence process in this 
system can be modelled using the composition of a phase- 
damping channel with a generalized amplitude damping 
channel: 

PAB{t) ^ £pO SaipAs), 

where the quantum operations £„ {n = P, 1) are written 
in the operator-sum representation as shown in Section 
1, i.e. 



t 



where K"{t) are the well-known Kraus' operators. For 
the generalized amplitude- damping channel the Kraus' 
operators are given by 



^0 




where, in the NMR context, 



and 



1-e 



p = l-exp(-Vri), 



(37) 



(38) 



with Ti being the longitudinal relaxation time of the 
qubit under consideration. In our case, each qubit have 
distinct Larmor frequencies, leading to also distinct relax- 
ation times. The measured spin-lattice relaxation times 
are 

TiCE) = 2.5 s and Ti("C) = 7 s. 

The Kraus' operators for the phase-damping channel are 
given by Eqs. ([Ml) with p = 1 — exp {—t/T2), where the 
measured transverse relaxation times associated with the 
the two qubits are given by 

T2(iH) = 1.8 s and T2("C) = 0.29 s. 

As no refocusing pulse was used, the effective transverse 
relaxation times are 

T^Cn) = 0.31 s and T^C^C) ^ 0.12 s. 

To verify the sudden-change in the decay rate and the 
robustness of correlations to phase-noise environments, a 
suitable initial state was prepared by mapping the NMR 
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FIG. 11. Description of pulse sequence used to implement 
the dynamics of quantum and classical correlations under de- 
coherence. The bottom (top) sequence corresponds to the 
representation of pulses applied to the (^^C) nuclei. At 
first, the SMP technique and a z-gradient pulse are applied 
to implement a diagonal state. In the sequence, a pseudo- 
EPR gate is performed. In the third step, the nuclear system 
and environment interact resulting in decoherence of quan- 
tum states. Finally, the read-out of the quantum state is per- 
formed for time intervals of m/4J (m — 0, 1, 2, 250) using 
the quantum state tomography procedure. 



deviation matrix Ap into a Bell-diagonal state (Eq. (HH)) 
with |ci|, |c2| > Icsj. An illustration of the pulse sequence 
utilized for accompanying the correlations dynamics is 
shown in Fig. [TlJ Using the deviation matrices obtained 
from quantum state tomography, the correlations were 
numerically computed for each time step to/4 J. The re- 
sults are presented in Fig. [T21 where the sudden-change 
in the decay rate of correlations is clearly present even 
in the presence of the strong thermal noise produced by 
the fluctuating fields of the nuclear spin environments. 
As already discussed the classical correlation should be 
constant under a phase damping, but a small decay is ob- 
served in Fig. [T^lwhich is due to the presence of the ther- 
mal environment. We see that the theoretical predictions 
presented are in good agreement with the experiment. 
We also verify the transition between two decoherence 
regimes [3l|. In the first regime the decoherence affects 
more strongly the classical aspect of correlation, and the 
quantum correlations have a small decay rate. After the 
sudden-change point, the classical correlation becomes 
more robust against decoherence and the noise affects 
more the quantum discord. It is worthwhile to note that 
the experimental deviation matrices have small coher- 
ences, which oscillate during the evolution. This leads 
to the small oscillations in the spectral line intensities 
used to perform the state tomography and, because the 
tomography process involves solving coupled equations, 
this oscillatory behaviour causes the oscillations observed 
in the correlations. 



2. Experimental Implementation of a Classicality Witness 

As we discussed so far, the quantification of quantum 
correlations involves quantum state tomography. How- 




0.25 



FIG. 12. Sudden change in behaviour and robustness of clas- 
sical and quantum correlations under decoherence. The green 
downward triangles are the experimental data for the quan- 
tum mutual information, while the blue circles and orange 
upward triangles represent the classical and quantum correla- 
tions, respectively. The black lines are the theoretical predic- 
tions. The initial state is analogous to the state in Eq. (|19|) 
with |ci|, jc2| > jcsl. The correlations are displayed in units 
of (eVln2)bit. 



ever, in some situations it is enough to identify the na- 
ture of correlations in a system. To meet this goal one 
usually uses correlations witnesses. It turns out, as no- 
ticed in Section 2, that the space of classical states is 
not convex. For that reason, linear witnesses cannot be 
used to identify the quantumness of correlations in sep- 
arable states. In what follows we present an experimen- 
tal implementation of the nonlinear classicality witness 
introduced in Section III El . The experiment was per- 
formed using the same NMR apparatus utilized to verify 
the sudden-change phenomenon of quantum discord un- 
der decoherence [S^- The pulse sequence used to imple- 
ment the protocol for the classicality witness is sketched 
in Fig. [131 Three states were prepared, a thermal equi- 
librium state pt, a quantum correlated state pQc and a 
classically correlated stated pec- The measured values of 
the witness for three different states is presented in Table 



rVl The thermal state has Wo^ 



0.05, so this is assumed 



as the classicality cutoff limit of our experiments . For 
the classical state Wp^^ = 0.04, which is within the clas- 
sicality cutoff limit. The quantum-correlated state has 
WpQc = 3.13, far above the classicality bound 0.05. As a 
matter of fact, the classicality witness perfectly sorts out 
quantum and classically correlated states. Table ITVl also 
shows the quantum discord computed from the exper- 
imentally reconstructed deviation matrices, which is in 
agreement with the classicality witness. However quan- 
tum discord is obtained after quantum state tomography 
followed by numerical extremization procedures. 

At last, we investigated the decoherent dynamics of 
the witness. We let the system, prepared in the initial 
state pqci evolve freely during a time period t„ = n5t 
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FIG. 13. Description of pulse sequence used to implement 
the Classicality Witness. The bottom (top) sequence corre- 
sponds to the representation of the pulses applied to the 
(^^C) nucleus. At first, the SMP technique and a 2-gradient 
pulse are used to implement a diagonal state. In the second 
step a pseudo-EPR gate is implemented. The third step is 
devoted to the interaction of the system with the environ- 
ment to produce de decoherence of quantum states. Finally 
the read-out of the quantum state was performed at m/4J 
time intervals with m = 0, 1, 2, 11, using the tomography 
procedure. 



TABLE IV. Witness, quantum discord, and classical correla- 
tion measured in three different initial states: quantum corre- 
lated pQc, classically correlated pec and thermal equilibrium 
pT- The witness was measured directly by performing the se- 
quence of pulses depicted in Fig. 1131 while the classical cor- 
relation and the symmetric quantum discord was computed 
after full QST and numerical extremization procedures. The 
correlations are displayed in units of (e^/ln2)bit. 





PQC 


pec 


PT 


Witness 


3.13 


0.04 


0.05 


Quantum Discord 


4.02 


0.00 


0.00 


Classical Correlation 


2.09 


7.15 


0.00 



{6t = 55.7nis, n = 0, 1, 2, • • • , 11) under the action of the 
NMR environment. The experimental results are shown 
in Fig. [TH In this figure we also show the values of the 
correlations. We obtain a fairly good agreement between 
the witness values and the correlation quantifiers when 
identifying classical states. 

IV. CONCLUDING REMARKS 

Quantum information processing has the potential of 
providing great advances, not only on the foundations of 
physics, specially to quantum mechanics and thermody- 
namics, but also on the technological ground. Therefore, 
studying the basis of quantum information is of funda- 
mental importance. 

Quantum correlations is thought to be responsible for 
the gain offered by quantum mechanics when compared 
with classical protocols. Our understanding about what 
is quantum and what is classical in a quantum correlated 
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system have changed dramatically since the seminal EPR 
paper in 1935 f]\. Initially, we thought that quantum cor- 
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FIG. 14. Measured witness and computed correlations for 
PQC relaxed during a time interval t„ — nSt {5t = 55.7ms, 
n — 0,1,2,-- - ,11), before performing the witness measure- 
ment protocol. The red circles represent the witness, the 
downward gray triangles are the amount of classical correla- 
tion, and the blue upward triangles represents the symmetric 
quantum discord. The correlations are displayed in units of 
(eVln2)bit. 

relations were the same as quantum non-locality. After 
that, entanglement was raised to the status of quantum 
correlations. In the last decade we realized that there 
are quantum correlations present even in separable (non- 
entangled) states and that such correlations could be 
used to over-perform classical systems in quantum infor- 
mation protocols. The most popular quantifier for such 
kind of correlation is called quantum discord • 

In this review we discussed the theoretical and the ex- 
perimental aspects of quantum discord in nuclear mag- 
netic resonance systems at room temperature. After dis- 
cussing the dynamics of quantum discord under deco- 
herence, we presented a witness for such a correlation. 
The experimental investigations of both issues was then 
treated, with a brief review of control techniques in NMR 
as well as of quantum state interrogation, a procedure 
known as quantum state tomography. 
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